Recovering coherence from decoherence: a method of quantum state reconstruction 
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We present a feasible scheme for reconstructing the quantum state of a field prepared inside a 
lossy cavity. Quantum coherences are normally destroyed by dissipation, but we show that at zero 
temperature we are able to retrieve enough information about the initial state, making possible to 
recover its Wigner function as well as other quasiprobabilities. We provide a numerical simulation 
of a Schrodinger cat state reconstruction. 
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I. INTRODUCTION 



^ ' The reconstruction of quantum states is a central topic in quantum optics and related fields j^,^ . During the past 
\ years, several techniques have been developed, for instance, the direct sampling of the density matrix of a signal 
mode in multiport optical homodyne tomography j^, tomographic reconstruction by unbalanced homodyning [Q, 
cascaded homodyning Q and reconstruction via photocounting There are also proposals of measurement of the 
electromagnetic field inside a cavity 0] as well as the vibrational state of an ion in a trap |^ . The full reconstruction of 
nonclassical field states as well as of (motional) states of an ion have been already experimentally accomplished. 
^ , The quantum state reconstruction is normally achieved through a finite set of either field homodyne measurements, or 
' selective measurement of atomic states Q in the case of cavities. This makes possible to construct a quasidistribution 
^0 (such as the the Wigner function) which constitutes an alternative representation of the quantum state of the field. 
' Nevertheless, in real experiments, the presence of noise and dissipation has normally destructive effects. In fact, 
\ as it has been already pointed out, the reconstruction schemes themselves also indicate loss of coherence in quantum 
^\ i: systems pO[ |. regarding this subject, a scheme for compensation of losses in quantum-state measurements has been 
0^ ' already proposed [pl|, and the relation between losses and s-parameterized quasiprobability distributions has been 
, already pointed out in [Q. The scheme on loss compensation in applies to photodetector losses, and consists 
^ ^ essentially of a mathematical inversion formula expressing the initial density matrix in terms of the decayed one. 
^ I ^ ' Our scheme, as discussed in fl^ , involves a physical process that actually enables us to store information about all 
r-| the quantum coherences of the initial state in the diagonal elements (photon distribuition) of the density matrix 
of a transformed state. By storing this information in the diagonal elements, it becomes much more robust under 
^ , dissipation, allowing us to recover the Wigner function of the initial field state in a time scale of the order of the 
energy decay time that is, of course, much longer than the extremely fast decoherence time scale that is normally 
associated with the dissipation of quantum coherences. 

We consider a single mode high-Q cavity where we suppose that a (nonclassical) field state /5(0) is previously 
^ . , prepared. The first step of our method consists in driving the generated quantum state by a coherent pulse. The 
' reconstruction of the field state may be accomplished after turning-off the driving field, i.e., at a time in which the 
cavity field has already suffered decay. We use the fact that by displacing the initial state (even while it is decaying) 
we make its quantum coherences robust enough to allow its experimental determination, at a later time, despite 
dissipation. We then show that the evolution of the cavity field is such that it directly yields the Wigner function 
of the initial nonclassical field simply by measuring the photon number distribution of the displaced field. For that 



we make direct use of the series representation of quasiprobability distributions |14|. A numerical simulation of our 
method is presented, and we take into account the action of dissipation while driving the initial field. 

This manuscript is organized as follows: in Sec. II we discuss, taking into account losses, the process of displacement 
of the initial field. In Sec. Ill we show how to reconstruct the initial cavity field after allowing the displaced field to 
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decay. In Sec. IV we present a simulation of the reconstruction of a Schrodinger cat state. In Sec. V we summarize 
our conclusions. 



II. DRIVING THE INITIAL FIELD 



We assume that the initial nonclassical field /5(0) is prepared inside a high Q cavity. The master equation in the 
interaction picture for the reduced density operator p relative to a driven cavity mode, taking into account cavity 
losses at zero temperature and under the Born-Markov approximation is given by |l5[| 



^ = - ^ [i/d, /5] + I {2apa) -a^ap- pa) a) 



with 



(1) 



(2) 



where a and are the annihilation and creation operators, 7 the (cavity) decay constant and a the amplitude of the 
driving field. 

We define the superoperators iZ and jC by their action on the density operator m& 



and 



It is not difficult to show that 



iZp = {a*a — aa^)/5 — p{a*a — aa^), 



£p — "fapa) — {a) ap + pa) a) 



and the formal solution of Eq. (|T|) can then be written as [0 

p{t) ~ exp (7^ + C)t /5(0) = exp(£f) exp 
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p(0). 



After driving the initial field during a time t^, the resulting field density operator will read 

p{U) = e^'-psm^ 

where 

P(m = b\m^)m). 

and with 

1 _ f.itdl'i 



P = 2a- 
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(3) 
(4) 

(5) 
(6) 

(7) 
(8) 

(9) 



This means that if we drive the initial field while it decays, during a time t^, this is equivalent to having the field 
driven by a coherent field with an effective amplitude /3 given in Eq. (ph. 



III. THE RECONSTRUCTION METHOD 



The driving of the initial field is carried out during a time t^. This procedure will enable us to obtain information 
about all the elements of the initial density matrix from the diagonal elements of the time-evolved displaced density 
matrix only. As diagonal elements decay much slower than off-diagonal ones, information about the initial state stored 
this way becomes robust enough to withstand the decoherence process. We will now show how this robustness can be 
used to obtain the Wigner function of the initial state after it has already started to dissipate. Once the injection of 
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the coherent pulse is completed, the cavity field is left to decay, so that its dynamics will be governed by the master 
equation in Eq. (0) without the first (driving) term in its right-hand-side. Therefore, the cavity field density operator 
will be, at a time t, given by 

P/3(t)=e(^+^)*pM0), (10) 

with 

Jp = ^apa\ Lp = — i^a) dp + pa) a) . (11) 
The next step is to calculate the diagonal matrix elements of Pi3{t) — exp (J + L)t pp in the number state basis, 

or 

Hm\m) = ^^^E?" rn ("l^/^WI")' (12) 

where q = 1 — e^^* . 

Now we multiply those matrix elements by powers of the function 

X{t) = 1 - 2e^*. (13) 

If we sum the resulting expression over m, we obtain the following simple sum 

Fw^-Yl = - E(-l)"("l^^(/3)/5(0)^(/3)l^)- (14) 

TT TT 

m— n— 

The expression in Eq. (|lj) is exactly the Wigner function corresponding to p (the initial field state) |jl^ at the point 
specified by the complex amplitude /3. Therefore we need simply to measure the diagonal elements of the dissipated 
displaced cavity field Pm{P',t) = {m\pi3{t)\m) for a range of /3's, the transformation in Eq. (|l^ ) in order to obtain 
the Wigner function of the initial state for this range. We note that after performing the sum, the time-dependence 
cancels out completely, leaving us a constant Wigner function, as it should be. Therefore the initial state may be 
reconstructed, at least in principle, at an arbitrary later time. In practice, however, the decay of the field energy will 
impose a limitation on the times during which we will be able to measure the PmS. 

The next step in our scheme is to measure the field photon number distribution Pm{P',t). In a cavity, particularly 
in the microwave regime where there are no photodetectors available, experimentalists have been forced to use atoms 
instead to probe the intra-cavity field. One way of determining Pm is by injecting atoms into the cavity and measuring 
their population inversion as they exit after an interaction time t much shorter than the cavity decay time. It is 
convenient in this case to use three-level atoms in a cascade configuration with the upper and the lower level having 



the same parity and satisfying the two-photon resonance condition. The population inversion in this case is |18 

(n+l)(n + 2) 



Wia;t + T) = E^"(";^) 



(52 S'^ 



■ cos (2(5„At) 



(15) 



n=0 

where r„ = [A + x{n + 1)] /2, S'^ — -|-A^(n-|-l)(n-t-2), A is the atom-field detuning, x is the Stark shift coefficient, 
and A is the atom-field coupling constant. Now we take A = (two-photon resonance condition) and x = 0. It is a 

1/2 

good approximation to make [{n + 1)(« + 2)] « n + 3/2, so that the population inversion reduces to 

oo 

W{a; t + T)^Yl ^»("; cos ([2^ + 3] Ar) . (16) 

This represents the atomic response to the displaced field. In order to obtain P„j from a family of measured population 
inversions, we need to invert the Fourier series in Eq. (|l^), or 

2A f'^™"^' 

Pm{l3;t)^— dTT^(i + t)cos([2to + 3]At) . (17) 



We need a maximum interaction time Tmax = tt/A much shorter than the cavity decay time, and this condition implies 
that we must be in the strong-coupling regime, i.e. A ^ 7. Of course after some time the atoms sent through the 
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cavity will not be able to respond anymore to the decaying field, hampering our reconstruction scheme. Nevertheless 
this happens at a time scale of td ~ 1/A, which is normally much longer than a typical decoherence time. 

Our scheme is easily extended to other (s-parametrized [^) quasi-probability distributions, which also may be 
expressed as a series Q, 

9 °° / -I- 1 \ " 

ms) = --^Y.[^,) hp.I"). (18) 

For this we have to multiply the photon number distribution of the displaced (and dissipated) field by a generalization 
of the function in Eq. (p^), or 

X(s;i) = l + -. (19) 

s — 1 

This increases our possibilities of measuring quantum states. We have available several quasiprobabilities, and we 
may choose a more convenient function depending on the particular conditions of a reconstruction experiment. 



IV. RECONSTRUCTION OF A SCHRODINGER CAT STATE 



We show now how our method may be applied to a specific case, e.g., to a reconstruction of a Schrodinger cat 
state, represented by a quantum superposition of two coherent states having distinct amplitudes, \a) and \ — a). The 
density operator corresponding to such state is 

p(0) = [\a){a\ + I - a){-a\ + e^'^|a)(-a| + e-'^\ - a){a\] , (20) 

where (f) is a relative phase and M a normalization constant. Schrodinger cat states are very fragile under dissipation 
pl| , and therefore are specially suitable for exemplifying our scheme. 

In a real experiment, just after having prepared the field in the state in Eq. ( pO| ) by one of the conventional methods, 
(see reference ||2^ for instance), the cavity will be driven by a coherent field, say, of amplitude ao. After that, at a 
time to, N conveniently prepared three- level atoms should cross the cavity, so that we are able to assign a particular 
population inversion Wao(io) = Pe{to) — Pg{to) to that time. These sort of measurements for the same driving field 
amplitude ao will be repeated for a range of times to, ti, t2 ... t = to + ""/-^j so that we obtain a set of values for the 
population inversion in that time interval. Now through numerical integration (see Eq. ([l7|)), we obtain the required 
photon number distribution Pm(Q!o,t). This is exactly what we need in our reconstruction scheme. The next step is 
to multiply Pm by the terms x(0™ — (1 ^ 2e'''*)™ and sum over m. This directly yields, according to our scheme, the 
value of the Wigner function of the initial field in Eq. ( |20| ) at the point ao = xo + iyo in phase space. We remark 
that the convergence of the series in Eq. (|lj) is guaranteed, because the photon number distributions of physical 
states normally decreases very quickly as m increases, and the statistical erros also become small simply because we 
are dealing with diagonal elements of the density operator in the number state basis . We need then to repeat 
this procedure M times, for different values of the driving field amplitude ai, a2 ... aM, to be able to cover enough 
points in phase-space and obtain the whole Wigner function of the original field. 

We have produced a numerical simulation of the above mentioned steps. Our simulation is going to be illustrated 



by the Wigner functions themselves. At < = 0.0, a Schrodinger cat state (Eq. 20) is generated within the cavity. Its 
corresponding Wigner function is shown in Fig. 1, where we note the characteristic interference structure. After the 
field has decayed, say, at t = O.I/7, the loss of coherence is indicated by the significant reduction of the interference 
structure, as seen in Fig. 2. Dissipation brings out the initially pure state very close to a statistical mixture 
We may instead drive the cavity with a coherent pulse (duration t^), at t = 0, in order to start our reconstruction 
procedure. After following the steps described above, we obtain the reconstructed Wigner function, as shown in Fig. 
3, which is essentially the one in Fig. 1. We note that both peaks as well as the interference structure characteristic 
of a Schrodinger cat state are entirely preserved. In our simulation we have considered fluctuations that might be 
present during the measurement of the atomic inversion. There might be experimental errors from various sources, 
such as flucutuations in the amplitude of the driving field as well as in the generated nonclassical state, which would 
cause distortions in the atomic inversion. In Fig. 4 we show the atomic inversion as a function of time and for a 
given value of the driving field amplitude (3 = (0, 2). We have obtained the field's Wigner function shown in Fig. 3 
departing from a family of those "distorted" atomic inversions, for different values of (3. 



4 



V. CONCLUSION 



In conclusion, we have presented a method for reconstructing the Wigner function of an initial nonclassical state 
at times when the field would have normally lost its quantum coherence. In particular, even at times such that the 
Wigner function would have lost its negativity, reflecting the decoherence process. A crucial step in our approach 
is the driving of the initial field immediately after preparation, which stores quantum coherences in the diagonal 
elements of the time evolved displaced density matrix, making them robust. We have therefore shown that the initial 
displacement transfers the robustness of a coherent state pO| against dissipation to any initial state, allowing the full 
reconstruction of the field state under less than ideal conditions. 

A natural application of our method would be the measurement of quantum states in cavities, where dissipation 
is difficult to avoid. Moreover, the application of the driving pulse at different times after the generation of a field 
state, would allow the "snapshooting" of the Wigner function as the state is dissipated. This means that valuable 
information about the (mixed) quantum state as well as about the decay mechanism itself could be retrieved while 
it suffers decay. The possibility of reconstructing quantum states even in the presence of dissipation may be also 
relevant for applications in quantum computing. Loss of coherence associated to dissipation is likely to occur in those 
devices, and our method could be used, for instance, as a scheme to refresh the state of a quantum computer [ p3| in 
order to minimize the destructive action of the environment. 
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FIG. 1. Wigner function of a Schrodinger cat state with a = 2 and = at t = 0. 




FIG. 2. Wigner function of a Schrodinger cat state with a = 2 and ^ = at i = O.I/7. 




FIG. 3. Wigner function of a Schrodinger cat state reconstructed &t t = O.I/7, with a = 2 and ^ = 0. The Wigner function 
at t = has been recovered (see Fig. 1). 
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FIG. 4. "Measured" atomic inversion (dashed line) and "exact" atomic inversion of (continuous line) the displaced cavity 
field as a function of time for a displacement amplitude /3 = (0, 2). 
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